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Abstract 

We study the inclusive decays of Df mesons, using data collected near the D* + Dj peak pro- 
duction energy E cm = 4170 MeV by the CLEO-c detector. We report the inclusive yields of Df 
decays to K+X, K~X, K° S X, ir+X, tt~X, tt°X, r/X, rj'X, <f>X, coX and / (980)X, and also decays 
into pairs of kaons, Df — > KKX. Using these measurements, we obtain an overview of Df decays. 
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The Df meson, consisting of a c and s quark, is the least extensively studied of the 
ground state charmed mesons. Here we present measurements of many inclusive yields from 
Df decay, thereby obtaining an overview of decays. 

Studies of inclusive branching fractions provide strong constraints on Monte Carlo sim- 
ulation. On completion of the measurements described here, we retuned our Monte Carlo 
decay table. The comparisons of Monte Carlo and data yields and spectra given below are 
after this retuning. 

In addition to providing an improved Monte Carlo decay table, our results allow some 
comparisons with expectations. 

Data for this analysis were taken at the Cornell Electron Storage Ring (CESR) using the 
CLEO-c general-purpose solenoidal detector, which is described in detail elsewhere [3,11, H, 
5J. The charged particle tracking system covers a solid angle of 93% of Att and consists of a 
small-radius, six-layer, low-mass, stereo wire drift chamber, concentric with, and surrounded 
by, a 47-layer cylindrical central drift chamber. The chambers operate in a 1.0 T magnetic 
field and achieve a momentum resolution of ~0.6% at p =1 GeV/c. Photons are detected 
in an electromagnetic calorimeter consisting of 7800 cesium iodide crystals and covering 
95% of in, which achieves a photon energy resolution of 2.2% at E 1 =1 GeV and 6% at 
100 MeV. We utilize two particle identification (PID) devices to separate charged kaons 
from pions: the central drift chamber, which provides measurements of ionization energy 
loss (dE/dx), and, surrounding this drift chamber, a cylindrical ring-imaging Cherenkov 
(RICH) detector, whose active solid angle is 80% of 4ir. The combined PID system has 
a pion or kaon efficiency > 85% and a probability of pions faking kaons (or vice versa) 
< 5% [5]. The detector response is modeled with a detailed GEANT-based p] Monte Carlo 
(MC) simulation, with initial particle trajectories generated by EvtGen 0] and final state 
radiation produced by PHOTOS [8]. The initial-state radiation is modeled using cross 
sections for D^DJ production at lower energies obtained from the CLEO-c energy scan [9| 
near the CM energy where we collect the sample. 

We use 586 pb _1 of data produced in e + e~ collisions at CESR near the center-of-mass 
energy s/s = 4170 MeV. Here the cross-section for the channel of interest, D* + D~ or Dtp*~, 
is ~1 nb [9j]. We select events in which the D* decays to D s ~f (94% branching fraction 10]). 
Other charm production totals ~7 nb |9j, and the underlying light-quark "continuum" is 
about 12 nb. 

Here we employ a double-tagging technique. Single-tag (ST) events are selected by fully 
reconstructing a D~ , which we call a tag, in one of the following three two-body hadronic 
decay modes: D~ — > KgK~, D~ — > 07r~ and D~ — > K*°K~. (Mention of a specific mode 
implies the use of the charge conjugate mode as well throughout this paper.) Details on 
the tagging selection procedure are given in Ref. 111. The tagged D s candidate can be 



either the primary D s or the secondary D s from the decay D* — > 7.D S . We require the 
following intermediate states to satisfy these mass windows around the nominal mass fioj ]: 
K° s -> 7T+7T- (±12 MeV), -> K+K- (±10 MeV) and K*° -> K+n' (±75 MeV). All 
charged particles utilized in tags must have momenta above 100 MeV/c to eliminate the soft 
pions from D*D* decays (through D* — > ttD). 

We use the reconstructed invariant mass of the D s candidate, M(D S ), and the mass re- 



coiling against the D s candidate, M recoi i(D s ) = y (E — E Ds ) 2 — (p — PdJ 2 , as our primary 
kinematic variables to select a D s candidate. Here (Eq, po) is the net four-momentum of the 
e + e~ beams, taking the finite beam crossing angle into account, pe> 3 is the momentum of the 



D s candidate, Ed s = Jm 2 D + p 2 D , and mr> s is the known D s mass [lOj. We require the re- 
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FIG. 1: The mass difference AM(D S ) = M{D S ) — rriD a distributions in each tag mode. We fit 
the AM(D S ) distribution (open circle) to the sum (solid curve) of signal (double-Gaussian) plus 
background (second degree polynomial, dashed curve) functions. 



coil mass to be within 55 MeV of the D* mass [10] • This loose window allows both primary 
and secondary D s tags to be selected. We also require a photon consistent with coming 
from D* s — > decay, by looking at the mass recoiling against the D s candidate plus 7 

system, M recoi i(D s 7) = J (E~ — E Da — E^) 2 — (p — p Ds — p 7 ) 2 . For correct combinations, 
this recoil mass peaks at mo s , regardless of whether the candidate is due to a primary or a 
secondary D s . We require |M recoil (Z) s 7) — m Ds \ < 30 MeV. 

The invariant mass distributions of D s tag candidates for each tag mode are shown Fig. [TJ 
We use the ST invariant mass sidebands to estimate the background in our signal yields 
from combinatorial background under the ST mass peaks. The signal region is \AM(D S )\ < 
20 MeV, while the sideband region is 35 MeV < |AAf(D a )| < 55 MeV, where AM(D S ) = 
M(D S ) — nio s is the difference between the tag mass and the nominal mass. To find the 
sideband scaling factor, the AM(D S ) distributions are fit to the sum of double-Gaussian 
signal plus second-degree polynomial background functions. We have 18586 ± 163 ST events 
that we use for further analysis. 

In each event where a tag is identified, we search for our signal inclusive modes recoiling 
against the tag. Charged tracks utilized in signal candidates are required to satisfy criteria 
based on the track fit quality, have momenta above 50 MeV/c, and angles with respect to the 
beam line, 9, satisfying | cos#| < 0.80. They must also be consistent with coming from the 
interaction point in three dimensions. Pion and kaon candidates are required to have dE/dx 
measurements within three standard deviations (3a) of the expected value. For tracks with 
momenta greater than 700 MeV/c, RICH information, if available, is combined with dE/dx. 
Candidate positrons (and electrons), selected with criteria described in Ref. 13], are required 
to have momenta of at least 200 MeV/c. 

For D+ -> K + X, D+ -> K~X, D+ -> vr+X, and D+ -> n~X modes, we count the 
numbers of charged kaons and pions recoiling against the tag where the tags are selected 
from both M(D S ) signal and sideband regions. Thus the combinatoric background is sub- 
tracted by using M(D S ) sideband events. The particle misidentification backgrounds among 
e, 7r and K are estimated by using the momentum- dependent particle misidentification rates 
determined from Monte Carlo and the e, 7r and K yields. Our identification can not distin- 
guish between muons and pions. So, we assume the muon yield equals the electron yield, 



and subtract accordingly. For D~ s 
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FIG. 2: Charged and neutral kaon and pion momentum spectra after background subtractions and 
efficiency corrections: (a) L>+ -> K + X, (b) Df -> vr+X, (c) D+ -► if" A", (d) D+ -► tt^X, (e) 
— > i^SX, (f) — > tt°X. The points are obtained from data and solid line indicates the Monte 
Carlo after tuning. Good agreement between data and tuned Monte Carlo is found. Monte Carlo 
is normalized to data based on tag yield. 



decay as a background and subtract it based on K s yields. The momentum-dependent (50 
MeV bins) efficiencies for track finding, track selection criteria, and particle identification 
are obtained from Monte Carlo simulation. 

The K s candidates are reconstructed in K s — ► 7r + 7r~ decay. The two pions have no PID 
requirements, and a vertex fit is done to allow for the K s flight distance. We identify 7r° 
candidates via tt° — > 77, detecting the photons in the Csl calorimeter. We require that 
the calorimeter clusters have a measured energy above 30 MeV, have a lateral distribution 
consistent with that from photons, and not be matched to any charged track. The Kg 
(or 7T°) yield is extracted by defining a signal region and sideband regions in the invariant 
mass distribution of the pion (or photon) pair. The sideband scaling factor is obtained from 
Monte Carlo, thus allowing for a non-linear background shape. We treat 7r° from K s decay 
as a background for the decay Z?+ — > n X and subtract it based on K s yields. 

The momentum spectra after all background subtractions and efficiency corrections are 
shown in Fig. [2] 

For the rj we use the 77 final state, which has a large branching fraction in 77 decays. To 
better handle the mild dependence of efficiency on rj momentum, we separate the 77 sample 
into two momentum ranges to measure the inclusive yields, one below 300 MeV/c and the 
other above. The 77 signal and background yields are determined by fits to a Crystal Ball 



function [13j, to account for the peak and the low mass tail, and background polynomial. 
We reconstruct rj' candidates in the the decay mode 77' — > 7r + 7r~77 with the 77 subsequently 
decaying into 77. Candidates for 77' are selected by combining 77 candidates within 3 r.m.s. 
widths of the nominal rj mass, with a pair of tt + 7t~. The mass difference between 7]7t + tt~ 
and 77 is then examined and fit to a Gaussian signal function and a background polynomial 
to extract the 7/ yields. The <p candidates are reconstructed in cf) — > K + K~ decay. We 
break the <fi sample into several momentum regions (200 MeV/c bins) since the <fi efficiency 
changes substantially with momentum. In each momentum region, the signals are fit with a 
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FIG. 3: Invariant mass distributions: (a) D+ -> rjX, (b) D+ -> rfX, (c) D+ -> 0X, (d) L>+ -> wl. 



sum of two Gausssian shapes and the background is fit to a polynomial. We reconstruct uo 
candidates in u — > 7r + 7r~7r° decay and extract the u signal yields from the 7r + 7r _ 7r° invariant 
mass distribution. The invariant mass distributions of 77, r)', 0, and cj candidates, summed 
over all momenta, are shown in Fig. [3j 

We form / (980) candidates using tt + tt~ pairs, / (980) — > tt + tt~. The pions are subject 
to the standard pion PID requirements. We find no significant evidence for the decay 
— > /o(980)X. We fit the invariant mass distribution of tt + tt~ pairs to a Gaussian signal 
function plus a second-degree polynomial background function and we obtain a yield of 30 ± 
47. The 90% confidence level upper limit is B(D+ -> / (980)X)B(/ (980) -> tt+tt") < 1.1% 
(statistical uncertainty only). Systematic errors are 6.8% for the efficiency estimation, 5.6% 
for the signal and background shape parameters, and other smaller errors, leading to a 
combined relative systematic error of 8.8%. We conservatively increase the upper limit by 
1.28 times the combined systematic errors, giving a upper limit, including systematic errors, 
of B(D+ -> / (980)X)/3(/o(980) -> tt+tt-) < 1.3%. 

We also measure the inclusive yields of _D + mesons into two kaons. After a tag is identi- 
fied, we search for the best kaon pair, based on particle identification likelihood or Kg mass, 
per mode recoiling against the tag. The kaon pair modes can be any of K S K S , K S K + , 
K S K~, K+K~, K+K+ or K~K~. For D+ -> K° S K+X and D+ -> K S K~X, we apply 
the sideband subtraction on K s candidate invariant mass distribution to remove the non- 
resonant decay background and get the signal yields. The _D + — ► K^K^X signal yield is 
extracted by defining a signal region on the scatter plot for the two Kg candidate invariant 
masses. In order to account for — > Kgir + 7C^X and — > tt + tt-tt + tt-X entering into 
the signal region of Df — > KgKgX, we perform a background subtraction which has two 
components. For all two charged kaons modes, we count the event numbers where at lease 
two charged kaons are found recoiling against the tag. In order to subtract the combinatoric 
background, we repeat the same procedure for each mode where the tags are selected from 
M(D S ) sidebands. The other possible backgrounds from generic D s decay are studied using 
Monte Carlo and found to be negligible. 

The double-tagging technique allows us to measure the inclusive yields for the decay 
— > K\X without directly detecting the K\. Instead, we reconstruct all particles in the 
event except the single K \ and infer the presence of a K\ from the missing four-momentum. 
Our signal is a peak in the missing mass squared distribution at the K\ mass squared. 
Similar missing-mass-squared techniques are used for —>■ K^KgX, — ► K < [K + X and 
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TABLE I: D s inclusive yield results. Uncertainties are statistical and systematic, respectively. The 
inclusive K® results are only used as a check for K s . The Df — > K\X yield requires a correction 



before comparing with the Df — > K S X yield, as explained in the text. PDG 10] averages are 
shown in the last column, when available. 



Mode 




Yield (%) 


K Q L Mode 


Yield(%) B{ 


n-J- 

Df ~ 


>7T + X 


119.3 ± 1.2 ± 0.7 








Df~ 




43.2 ± 0.9 ± 0.3 








Df~ 


>7T°X 


123.4 ± 3.8 ± 5.3 








Df~ 


>K+X 


28.9 ± 0.6 ± 0.3 






20 


Df~ 


>K~X 


18.7 ± 0.5 ± 0.2 






13 


Df~ 


>r,X 


29.9 ± 2.2 ± 1.7 








Df - 


>r]'X 


11.7 ± 1.7 ± 0.7 








■m — s 1 

Df - 


xpX 


15.7 ± 0.8 ± 0.6 








Df~ 


>ujX 


6.1 ± 1.4 ± 0.3 








Df~ 


>/o(980)X,/ (980) - 


-» vr+vr- < 1.3% (90% CL) 








Df~ 


>K%X 


19.0 ± 1.0 ± 0.4 


Df 


-^K\X 


15.6 ± 2.0 20 


Df~ 


>K%K%X 


1.7 ± 0.3 ± 0.1 


Df 


^K° L K° S X 


5.0 ± 1.0 


Df~ 


>K S K+X 


5.8 ± 0.5 ± 0.1 


Df 


->KlK+X 


5.2 ± 0.7 


Df~ 


>K° S K-X 


1.9 ± 0.4 ± 0.1 


Df 


^K° L K-X 


1.9 ± 0.3 


Df~ 


>K+K-X 


15.8 ± 0.6 ± 0.3 








Df~ 


>K+K+X 


< 0.26% (90% CL) 








Df~ 


>K-R-X 


< 0.06% (90% CL) 









18 
14 
14 
12 



± 14 



Df — > K\K~X modes by requiring there must be a Kg, K + or K~ recoiling against the 
tag. Note that if the D s decay contains two or more K^s, we do not find any K\. Due to 
the low statistics and large systematic uncertainties, we quote the inclusive K\ results only 
as a check for Kg. 

The inclusive yields are listed in Table HI For the Kg modes, the corresponding K\ modes 
are listed as a comparison. The value of the decay Df — > K\X is only for Df decaying into 
a single K\. So one should not directly compare the values of Df — > KgX and Df — > K\X 
in Table [H One can correct the single K\ inclusive yield by adding two times the inclusive 
yield of Df K° L K° L X (assuming B(Df -> K° L K° L X) = B(Df -> K S K S X)). All the K9 
modes are consistent with K s modes. In the last column of Table [U we show PDG 10 
averages, when available. 

We have considered several sources of systematic uncertainty. The uncertainty associated 
with the efficiency for finding a track is 0.3%; an additional 0.6% systematic uncertainty for 
each kaon track is added [B| . The relative systematic uncertainties for ir° and K s efficiencies 
are 4.2% and 1.8%, respectively. Uncertainties in the charged pion and kaon identification 
efficiencies are 0.3% per pion and 0.3% per kaon [Bj]. All efficiencies from Monte Carlo have 
been corrected to include several known small differences between data and Monte Carlo 
simulation. 

The quark-level diagrams contributing to Df decay are shown in Fig. |H We classify 
"quark-level final states" as ss (as would come from Fig. Ufa)), s (Fig. H(b)), sss (Fig. BJc)), 
ss (Fig. lU^d)), and "no strange quarks" (Fig. @Je) and Fig. IH(f)). The ss final state is 
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FIG. 4: The typical Feynman diagrams of decays: (a) Cabibbo-favored decay, (b) single- 
Cabibbo-suppressed decay, (c) single-Cabibbo-suppressed decay, (d) double-Cabibbo-suppressed 
decay, (e) short-range annihilation decay, (f) long-range annihilation decay. 



Cabibbo-favored. The s and sss final states are singly-Cabibbo-suppressed, the ss final 
state is doubly-Cabibbo-suppressed, and the "no strange quarks" final state arises from 
short-range (Fig. 0Je)) and long-range (Fig. \Mf)) annihilation diagrams (While Fig. H)^f) 
shows the ss annihilating into gluons, here we also include its rescattering into uu or dd). 

The ss final state can hadronize as KKX, but also as r)X, rj'X, or (j)X. The s final state 
will hadronize as KX. The sss final state in principle can hadronize as KKKX, but there 
will be limited phase space for this, so Ki]X, KrfX, KifiX are probably more likely. The 
ss final state will hadronize as KKX, but being doubly-Cabibbo-suppressed, can probably 
be ignored. 

We have performed a global fit to our measurements. For this, we have branching fractions 
B(XX). In particular, for ss quark-level final states, we write B(D S — > ss) = B(ss), 
B(D S -> ss -> r]X) =_ B(r)), B{D S ^ ss -> rfX) = B{rf), B{D S -> ss -> <pX) = B(</>), 
and B(D S -> ss -> Ki^X) = B(A7iQ. Thus B(ss) = £(77) + S(^) + + Note 
that B(D S —>■ ss —>■ i]X) is the branching fraction for primary production of r] (not from 7/ 
decay), from the quark- level state ss. The free parameters in our fit are B{rj), B(r]'), B(<p), 
and B(KK), which we adjust to obtain the best fit. 

For the s quark-level final state, we note that B(D S — > s) = B{s) w (Kd/Ksl 2 x B(ss). 
Thus, we do not adjust B(s) in the fit, but write B(s) — C\ X iKd/Ksl 2 x B(ss), where C\ is 
a phase space correction factor, probably a bit larger than 1.0. We take C% to be 1.25 ±0.25. 

We break the sss quark-level final state into 4 separate pieces, as we have done with the 
ss final state. Thus B(D S — > sss) = B(sss) is made up of B(D S — > sss — > r/sX) = B(rjs), 
B(D S -> sss -> r/'sX) = B(77's), i3(D s -> sss -> 0sX) = B((j)s), and i3(D s -> sss -> 
i^/sTsX) = B(KKs). Thus £(sss) = i3(r/s) + B(^'s) + £(0s) + B(KKs). We note that 
£>(sss) ~ |Kis/Kid| 2 x S(ss). So again, we do not adjust any of the pieces making up B(sss), 
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but rather write 

B( V s) = C 2 x \V W /V nd \ 2 x B( V ) (1) 



B( V 's) = C 2 x llVKdl 2 x B{rf) (2) 

B(<f)s)=C 2 x \V m /V ud \ 2 x (3) 

= C 2 x |K S /Kd| 2 x B(KK) (4) 

The quantity C 2 , like C*i, is a phase space correction factor, expected to be smaller than 
1.0. We take it to be 0.75 ± 0.25. Assuredly the true phase space correction factors would 
be different for 77, 77', 0, and KK. We neglect this in our fit, allowing for it as a systematic 
error. 

For the doubly- Cabibbo-suppressed decays, we estimate B(D S — > ss) = B(ss) = C3 x 
KKd/KsXVus/Vud)! 2 x B(ss). This term is down a factor of 400 from the dominant term, 
and has essentially no effect on our fit. We take C3 = 1.0 ± 1.0. 

Finally, there are annihilation diagrams. We write ^(Annihilation) = B{D+ — > y^v) + 
B(D+ — > t + u) + B(Df — > Other Annihilation). One of our goals in performing the global 
fit is to get an estimate of B(D+ — > Other Annihilation). In our fit, we use B(Dt t + u) = 
(5.62 ± 0.41 ± 0.16)% [Il|, and B(D+ -> fj, + u) = (0.565 ± 0.045 ± 0.017)% Q. 

It is possible for a D s decay to contain more than one of 77, rf, (ft, KK, e.g. r\r\, r)<f), 
etc. From energy conservation, one of an allowed pair must be 77. So, we include a yield 
£>(extra 77) to allow for this. We searched for — > r/r/X, — ► rjrj'X, and — > r)<fiX. 
We found no clear signals, obtaining a summed yield of (6.0 ± 3.9)%. In our global fit, we 
take £>(extra rj) to be 6.0%, and include the ±3.9% in the systematic error. 

Another source of 77 and 77' is the quark-level decay D s — > s (Fig. BJb)). Here, the 77 
or 77' will come not from their ss component, but from their uu and dd components. At 
quark level, the decay is D s — > udds, so making rj or 77' is natural. We assume that this 
diagram gives an 77 a fraction fi of the time, and an 77' a fraction f 2 of the time, where 
/1 + ji < 1- While one can make quark-level predictions of what to expect for j\ and f 2 , 
we take the conservative position of allowing them the full range, < fi + f 2 < 1, and take 
fi = /2 = 1/4, in the middle of the allowed range. 

For our global fit, we write 



2 _ ^Y^-(g(77)+g(T ? 8)+g(77 , ->T7X)x(g(7 ? / )+g(77 / 8)+/ 2 xg(s))+g(extra gVfAxgfg)) x 2 



, Y^gfo / )+gfa'5)+/ 2 xg(g)K 2 
^ Y^(g(0)+g(#)) ^ 2 ! 



^ Y g Jc-(g(jrg)+g(JrA'a)+g(^g J K')x(g(0)+g(^3))+g(as)) ^2 
( 



5y 



+ 

AT if 



8y* 



1 K 



( 5 ) 

Here Yj is the central value of a measurement, and is the error on that measurement. 
As 7/ decays to 77, and <fi decays to KK, our \ 2 needs the branching fractions for those 
decays, B(rj' -> r]X) and B{<j> -> tfA"). We take these from PDG [Tj. Better than words, 
Eq. (jSJ) gives the meaning of the various B(XX) parameters. Thus, the measured yield of rj, 
Y v , has contributions from primary production of 77 from the ss quark state (B(i])), primary 
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TABLE II: Results from the global fit. The central values of parameters are listed in second 
column. The errors: Si is statistical uncertainty, 62 is from phase space factor C\ = 1.25 db 0.25, 
83 is from phase space factor C2 = 0.75 ± 0.25, £4 is from f\ + f2 = 0.5 ± 0.5, and £5 is from the 
£(extra 77) = (6.0 ± 3.9)%. 



Parameter 




Value(%) 




Error(%) 














S3 




£5 


B(D S -^ss- 


+ 7/X) 


14.7 


2.9 


0.2 


0.2 


1.0 


3.7 


B(D S -> ss- 


+ T/X) 


10.3 


1.7 


0.2 


0.1 


1.0 


0.1 


B{D S -> ss- 




15.1 


1.0 


0.0 


0.2 


0.0 


0.0 


B{D S -> ss- 




25.4 


1.2 


0.3 


0.6 


0.1 


0.1 


B{D S -> as) 




65.6 


2.7 


0.7 


1.0 


1.8 


3.5 


£>(Other Annihilation) 


21.5 


2.8 


0.1 


0.3 


2.0 


3.9 



production of 77 from the sss quark state (B(rjs)), primary production of 77 from the s quark 
state (/1 x B(s)), production of r\ from decay of 77', the rf being from the ss quark state 
(B(r]') x B(r]' — > T]X)), or the rj' being from the sss quark state (B(r)'s) x #(7/ — ► r]X)), 
or from the s quark state (/2 x B(s) x /3(r/' — > r)X)), and finally of "extra 77's", rj that 
accompanies an rj, rj', or already recorded (£> (extra rj)). The measured yields for rj' and 0, 
while not as complicated, have some of the same features. Note that, as described earlier, 
our measured yield of di-kaons, Ykk, includes KK and KK and KK pairs. There is a 
subtlety in the last line of Eq. (jSJ). The decay D s — ► sss always makes at least one kaon, 
and when the decay is D s — > KKs, i.e., B(KKs), makes 2 more. Line 5, for the kaon yield, 
properly handles this. 

We minimize \ 2 by varying /3(r/), B(rj'), B(<p), and B{KK). All other B(XX) parameters 
are fixed as previously described. Further, we have the unitarity requirement B(ss)+B(sss) + 
B{s) + B(ss) + ^Annihilation) = 1.0. Our fit gives B(r]), B(r)'), B{<f>), B(KK), and hence 
B(ss), B(sss), B(s), and £>(ss). Unitarity then gives /3(Other Annihilation). Results are 
given in Table [Til 

We have five measurements, and four free parameters. So it would appear that there 
is one degree of freedom. However, the single kaon and di-kaon measurements are highly 
correlated, so we effectively have more like four measurements. This is reflected in the x 2 
of the fit, which is 0.03. We have also made a fit leaving the di-kaon term out, and a fit 
leaving the single kaon term out. These fits give essentially the same result as the nominal 
fit with both terms included. 

In interpreting the results in Table [IT] it should be recognized that the decay products 
of the true "other annihilation" diagrams will include some D s — > gluons — > ss events, thus 
being treated as part of B(ss) rather than "other annihilation". Also, the gluons will make 
uu, dd, which will sometimes make 77, 77', again being treated as a contribution to B(ss). 
Thus £> (Other Annihilation) should be viewed as a lower bound, B(rj), B(r)'), B(cj)), B(KK) 
as upper bounds, on contributions from the various diagrams in Fig. |H On the other hand, 
an overestimate of £>(extra 77) will give an overestimate of /3(Other Annihilation). 

We can obtain a conservative lower bound on £>(Other Annihilation) by setting fo = f 2 = 
and /3(extra 77) = 0. That gives /3(Other Annihilation) = 13.3 ±3.0%, i.e., > 9.5% at 90% 
C.L.. 



10 



We use our measurements of the total kaon yield and the total di-kaon yield to get a 
measurement of the singly- Cabibbo-suppressed rate. If there were no tri-kaon events, then 
(total kaon yield) minus 2 x (total di-kaon yield) would give (single kaon yield) which would 
include the s final state, and that fraction of the sss final state for which the ss component 
hadronized as 77, r/, or 0. Tri-kaon events complicate the situation. As mentioned earlier, 
in counting di-kaons, a given charge pairing (K + K + , K + K®, K + K~ etc.) is counted once. 
Thus KgK^KgX is counted as one di-kaon, while K + KgKgX is counted as two, K + K° S K~ X 
as three. For the total kaon yield, a tri-kaon event is counted as 3 kaons, In taking (total 
kaon yield) minus 2 x (total di-kaon yield) as a way of counting singly-Cabibbo-suppressed 
yield, the "right" answer for a tri-kaon event is +1, and what we actually obtain is +1, — 1, 
and —3, for the different tri-kaon events, on average —1 instead of +1. Thus, our proposed 
procedure will underestimate the singly-Cabibbo-suppressed rate. To the extent that the 
tri-kaon rate is small, the underestimate is small. We estimate and apply a correction. 

Our numbers are: total kaon yield is (85.6 ± 2.3)%, total di-kaon yield is (39.9 ± 1.8)%. 
The errors are highly correlated. Taking correlations into consideration, we find kaon — 2 x di- 
kaon is (5.8 ± 2.2)%. Taking B(sss)/B(ss) to be ~ 1/20, and B(sss -> tri-kaon) /B(sss) 
to be < B(KK)/B(ss) = 0.39, our correction factor for the presence of tri-kaon decays is 
< (65.6x^x0.39x2)%. Thus, the correction factor is < 2.6%. Taking it to be (1.3 ± 1.3)%, 
the measured branching fraction for D s — > single-Cabibbo-suppressed is (7.1 ± 2.2 ± 1.3)%. 
The expected branching fraction is (|Vus/Kd| 2 + | V c d/V cs \ 2 ) x B(ss) « ^ x B(ss). Taking 
B(ss) from Table HTl we see fine agreement between expectations and measurements. 

From our global fit, we can compute the minimum yields of 7r + , 7r~, and tt° for each 
category. For example, for the Cabibbo-favored decay £)+ — > ss — > r]X, with 14.7% yield, 
we compute the yields of 7r + , 7r~, and 7r° that come from a 14.7% 77 yield. To this we add 
14.7% 7r + yield, since that must be present to conserve charge. (This is an overestimate, 
because semileptonic decays have charge conserved via e + or consequently we perform a 
subtraction to allow for that.) For —>■ sss — > rjsX, with 0.6% yield, similarly we compute 
the yields of 7r + , 7r~, and tt° that come from a 0.6% r\ yield. Charge conservation might be 
achieved by a ir + , but also by a K + . Lacking any information on how much comes from 
7r + , how much from K + , we assume half from each. Our global fit gives a single number 
B(KK) = 25.4%, for the di-kaon yield. To determine the 7r + , it , and ir~ yields, we need 
yields for the separate di-kaon combinations, KgKg, K$K + , K$K~ , etc. For our calculation, 
we take the measured di-kaon yields from Table [U and normalize them so their sum equals 
B(KK). (Where we have only an upper limit, we use half of it for the "measurement"). 

The results of our computation are given in Table IHIl There one sees that the yields of 
7r + , 7r~, and 7T° should be larger than 96.2%, 20.5%, and 46.8%, respectively. The observed 
yields are indeed larger than these numbers. Thus, on average, 1/4 of the D s decays will 
contain an additional 7r + 7r~ pair, and 3/4 of the D s decays will contain an additional 7r° (or 
1/2 contain one additional 7r°, 1/8 contain two additional 7r°'s). 

For the 21.5% yield of D s — > Other Annihilation decays, we know nothing about the pion 
content other than that there will be one 7r + to conserve charge. One might reasonably 
expect that a substantial fraction of the 1/4 of the D s decays containing an additional 7r + 7r~ 
pair would be in the "Other Annihilation" decays. As for the additional tt° in 3/4 of the 
decays, that can appear any place, e.g., as converting a charge-conserving n + into a p + . 
They will probably appear disproportionally in the "Other Annihilation" decays, as these 
start (in our table) with fewer particles. 

The inclusive uj yield, D s — > ujX, of 6.1 ± 1.4%, is substantial. While u has an ss 
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TABLE III: The minimum yields of ir + , it~ , and 7r° for each category. We compute the yields 
of 7r + , 7r~, and 7T° that come from signal particles. In addition to that, we add charged pions to 
conserve charge. Semileptonic decays have charge conserved via e + or u + , consequently we perform 
a subtraction to allow for that. 



Charge 
Conservation 



Particle 
Decay 



Total 
Yields 



Mode 




B (%) 


+ 


ix- 


+ 


7T~ 







7T~ 




7T 


Df - riX 




14.7 


14.7 


0.0 


4.0 


4.0 


17.7 


18.7 


4.0 


17.7 


Df - nsX 




0.6 


0.3 


0.0 


0.2 


0.2 


0.7 


0.4 


0.2 


0.7 


D+ - r/X 




10.3 


10.3 


0.0 


9.7 


9.7 


12.7 


20.0 


9.7 


12.7 


Dt rfsX 




0.4 


0.2 


0.0 


0.4 


0.4 


0.5 


0.6 


0.4 


0.5 


Dt - <px 




15.1 


15.1 


0.0 


2.4 


2.4 


2.5 


17.5 


2.4 


2.5 


Dt -► 0sX 




0.6 


0.3 


0.0 


0.1 


0.1 


0.1 


0.4 


0.1 


0.1 


D+ -> Extra nX 




6.0 


0.0 


0.0 


1.6 


1.6 


7.2 


1.6 


1.6 


7.2 


Dt — > (no 77, 77') 




2.1 


1.0 


0.0 


0.0 


0.0 


0.0 


1.0 


0.0 


0.0 


£>+ ^sl.X^yj 




1.0 


0.5 


0.0 


0.3 


0.3 


1.2 


0.8 


0.3 


1.2 


D+^sX,X^r/ 




1.0 


0.5 


0.0 


1.0 


1.0 


1.3 


1.5 


1.0 


1.3 


Dt - Xgi^M* 


3.3 


3.3 


0.0 


0.0 


0.0 


0.0 


3.3 


0.0 


0.0 




)* 


11.4 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


D+^K° S K-(K° L K- 


)* 


3.7 


7.5 


0.0 


0.0 


0.0 


0.0 


7.5 


0.0 


0.0 


D+ -► K+K-(-(p)X 




7.9 


7.9 


0.0 


0.0 


0.0 


0.0 


7.9 


0.0 


0.0 


Dt -► K+K+X 




0.1 


0.0 


0.1 


0.0 


0.0 


0.0 


0.0 


0.1 


0.0 


Dt -»• K-R-X 




0.03 


0.1 


0.0 


0.0 


0.0 


0.0 


0.1 


0.0 


0.0 


Dt - 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 






10.7 


-10.7 


0.0 


0.0 


0.0 


0.0 


-10.7 


0.0 


0.0 


L>+ -► r+zy 




5.6 


0.0 


0.0 


4.1 


0.8 


2.9 


4.1 


0.8 


2.9 


£>+ - 




0.6 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


— ► Other Annihilation 


21.5 


21.5 


0.0 


0.0 


0.0 


0.0 


21.5 


0.0 


0.0 


Minimum Yields 
Observed Yields 
Additional Yields 
















96.2 

119.3 

23.0 


20.5 
43.2 
22.7 


46.8 

123.4 

76.7 



component, it is very small, so it is unlikely that very much of the uj yield comes from the 
ss component of Dt — > ssX. At quark level, this is Df — > ssud, and a decay Dt — > n + r]uj 
is quite possible. A decay Dt — > 7i + i]'uj, from energy considerations, is just barely possible. 
From the decay Dt — > sss, uj could come from Dt — > K + i]uj (barely), but not from 
Dt — > K + i]'uj. From Dt — > sX, it can come from Dt — > K + ujX, with lots of phase 
space. And from "Other Annihilation" , there are lots of possibilities. In summary, with the 
data we now have in hand, we can not say much about the origin of the 6% uj yield. A 
search for Dt exclusive decays will be reported in a separate paper. (We should note that 
our inclusive uj measurement came towards the end of the work described here, and so was 
not included in the retuning of the Monte Carlo decay table. CLEO's D s Monte Carlo decay 
table produces far fewer cj's than the 6% we observe.) 

In summary, we report several measurements of Dt inclusive decays with significantly 
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better precision than current world averages. 
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